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ABSTRACT: We study the scalar, vector and tensor two-gluon and trigluon glueball spectra 
in the framework of 5-dimension dynamical holographic QCD model, where the metric 
structure is deformed self-consistently by the dilaton field. For comparison, the glueball 
spectra are also calculated in the hard-wall and soft-wall holographic QCD models. In 
order to distinguish glueballs with even and odd parities, we introduce the positive and 
negative coupling between the dilaton field and glueballs, and for higher spin glueballs, we 
introduce a deformed 5-dimension mass. With this set-up, there is only one free parameter 
from the quadratic dilaton profile in the dynamical holographic QCD model, which is fixed 
by the scalar glueball spectra. It is found that the two-gluon glueball spectra produced in 
the dynamical holographic QCD model are in good agreement with lattice data. Among 


six trigluon glueballs, the produced masses for 157 and 277 are in good agreement with 
lattice data, and the produced masses for 0~~, 0*~ and 2*7 are around 1.5 GeV lighter 
than lattice results. This result might indicate that the three trigluon glueballs of 077, 
0*~ and 2*~ are dominated by three-gluon condensate contribution. 
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1 Introduction 


Quantum chromodynamics (QCD) is accepted as the fundamental theory of describing 
strong interaction. In the high energy regime, QCD has the property of asymptotic freedom 
and perturbative QCD calculations have been tested with high precision. However, in the 
low energy regime, the nonperturbative aspect related to QCD vacuum properties and 
hadron spectra remains as outstanding challenge. The nonabelian feature of QCD makes 
it possible to form bound states of gauge bosons, i.e. glueballs (gg, ggg, etc. ) [1]. The 
gauge field plays a more important dynamical role in glueballs than that in the standard 
hadrons, therefore study particles like glueballs offers a good opportunity of understanding 
nonperturbative aspects of QCD. 

The glueball spectrum has attracted much attention more than three decades [1], and 
it has been widely investigated by using various non-perturbative methods. For example, 
from first principles calculation by using lattice QCD [2-6], by using flux tube model [7] 
as well as by using QCD sum rules [8-10]. For more information, please refer to review 
papers [11]. 

The discovery of the gravity/gauge duality, or anti-de Sitter/conformal field theory 
(AdS/CFT) correspondence [12-14] offers a new possibility to tackle the difficulty of 
strongly coupled gauge theories, for reviews see Ref. [15]. In recent decades, many efforts 
have been invested from both top-down and bottum-up approaches in examining nonper- 
turbative QCD properties, e.g., QCD equation of state, phase transitions, fluid properties 
of quark-gluon plasma [16], meson spectra [17-19], baryon spectra [20], as well as in the 
glueball sector [21-25]. It is expected that the holography approach can shed some light 
on our understanding of the nonperturbative aspects of QCD. 

QCD is a non-conformal gauge theory, and the Sakai-Sugimoto (SS) model [26] is one 
of the most successful non-conformal top-down holographic QCD models. The glueball 
spectra in the Sakai-Sugimoto model has been investigated in literatures, see Ref. [27]. 
Glueballs have also been widely studied by using the bottom-up approach [23], where most 
studies are based on hard-wall [17] and soft-wall holographic QCD models [18] with the 
conformal AdS; background metric. 

A successful holographic QCD model should grasp two main features of nonpertur- 
bative QCD properties, i.e. the spontaneous chiral symmetry breaking and color charge 
confinement. The dynamical holographic QCD model (DhQCD) which can describe both 
chiral symmetry breaking and confinement has been constructed in Ref. [28-30]. In this 
model, the gluon dynamics background is determined by the coupling between the graviton 
and the dilaton field ®(z), which is responsible for the gluon condensate and confinement, 
and the scalar field X(z) is introduced to mimic chiral dynamics. Evolution of the dilaton 
field and scalar field in 5D resemble the renormalization group from ultraviolet (UV) to 
infrared (IR). This DhQCD model describes the scalar glueball spectra and the light meson 
spectral quite well [28-30]. Further studies [31-33] show that this DhQCD model can also 
describe QCD phase transition, equation of state of QCD matter and temperature depen- 
dent transport properties, including shear viscosity, bulk viscosity, electric conductivity as 
well as jet quenching parameter. 


For the scalar glueball spectra, it was shown in Ref.[29] that, comparing with the results 
in the hard-wall and soft-wall holographic QCD models [23], the scalar glueballs including 
the lowest state and excited states can be surprisingly well described in the DhQCD model 
. However, the scalar glueball 0** has the same quantum number as the scalar quarkonium 
qq and tetraquark gqqq [34] states, and the complexity of determining the glueball states lies 
in that gluonic bound states might always mix with gq and qqqq states. For example, one 
has to distinguish the lightest scalar glueball state among 19 scalar mesons observed in the 
energy range below 2 GeV [35, 36]. Therefore, it is interesting to investigate odd glueballs 
with unconventional quantum numbers which cannot be carried by quark-antiquark bound 
states. These include JPC = 0-~,0+~,2+~,3-* glueballs, which can only be made of 
at-least three-gluon bound states. 


This motivate us to investigate the whole glueball spectra including (scalar, vector as 
well as tensor glueballs and their excitations) in the framework of the DhQCD model. The 
paper is organized as following: In Sec.2 we give the operators of two-gluon and trigluon 
glueballs. We introduce the dynamical soft-wall holographic QCD model in Sec.3, and 
calculate the the glueball spectra in the dynamical holographic QCD model in Sec. 4, it 
is found that higher-spin glueballs are very heavy comparing with lattice data, and the 
even and odd parities cannot be distinguished. Therefore, we introduce a deformed 5- 
dimension mass for higher spin glueballs, and in order to distinguish glueballs with even 
and odd parities, we introduce the positive and negative coupling between the dilaton field 
and glueballs.With this set-up, we calculate the glueball spectra in the modified dynamical 
holographic QCD model in Sec. 5 and find that the two-gluon glueball spectra are in good 
agreement with lattice data and the trigluon glueball spectra agree with results from QCD 
sum rules. Finally, a short summary is given in Sec.6. 


2 Two-gluon and trigluon glueballs 


The AdS/CFT correspondence establishes a one-to-one correspondence between a certain 
class of 4D local operators in the VV = 4 superconformal gauge theory and 5D supergravity 
fields representing the holographic correspondents in the AdSs x $° bulk theory. According 
to AdS/CFT dictionary, the conformal dimension of a (f-form) operator on the ultraviolet 
(UV) boundary is related to the M? of its dual field in the bulk as follows [12-14] : 


M3 =(A- f\(A+f-4). (2.1) 


In the bottom-up approach, for example in the holographic QCD models, one can 
expect a more general correspondence, i.e. each operator O(a) in the 4D field theory 
corresponds to a field O(a, z) in the 5D bulk theory. To investigate the glueball spectra, 
we consider the lowest dimension operators with the corresponding quantum numbers and 
defined in the field theory living on the 4D boundary. We show the two-gluon and trigluon 
glueball operators and their corresponding 5D masses in Table 1. 
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Table 1. 5D mass square of two-gluon and trigluon glueballs.The operators are taken from [22] 
and [9, 10]. 


For trigluon glueball 077, the detailed structure of the operator is given in Ref. [9] 
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where d% stands for the totally symmetric SU,(3) structure constant and Top = Jap — 
ag /0?. 
The interpolating currents of the 2+7 oddball 2+~ takes the form as [10], 


Jaa (2) =93d""[G4,(2)] [Gb p(@)IGea(a)] (2.6) 
Sno?” (2)=93d""[Gt, (2)IIGrp(@)Gba(a)] ; (2.7) 
Sua“ O= a (Gt, (@)[Gep(#)][Gea(@)] (2.8) 
Sno?” (2)=95 d” (GE, (2) [Gp (2) [Ga(2)] - (2.9) 


3 The dynamical soft-wall holographic QCD model and gluodynamics 


The dynamical soft-wall holographic QCD model is described in Ref.[29]. The pure gluon 
part of QCD can be modelled by the 5D graviton-dilaton coupled action: 
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where Gs is the 5D Newton constant, gs, ® and V4 are the 5D metric, the dilaton field 
and dilaton potential in the string frame, respectively. The metric is chosen to be 


Jun = b2(z)(dz? pac dn’), bs(z) = eAs(), (3.2) 
Under the conformal transformation 


gun = 9iune 2°, VE = PVE, (3.3) 


Eq.(3.1) can be rewriten in the Einstein frame 
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The Einstein equations are 
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Substituting the metric of Eq.(3.2) into the above equations, we can obtain: 
A t9 4 19 
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bp(z) = b.(z)e~32 = e420), Ap(z) = A,(z) — mo (3.9) 
In the string frame, the above two equations of motion are 
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We take the same dilaton field as that in the KKSS model or soft-wall holographic 
QCD model [18], i.e., 


Oo = pre. (3.12) 

It is simple to solve the metric Ag and the dilaton potential VĚ (®) in the quadratic dilaton 
background 
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with oF (a; z) the hypergeometric function. 


4 Glueball spectra in the dynamical soft-wall holographic QCD model 


4.1 Scalar glueballs 


The 5D action for the scalar glueball Y(z, z) in the string frame takes the form as that in 
the original soft-wall model [24, 25} 


Sg = - | ëz Tage [OuGO"G + M3 59"). (4.1) 


It is notice that the metric structure in the dynamical soft-wall model is solved from Eq. 
(3.8) instead of AdSs. 
The Equation of motion for the scalar glueballs Y is given below as 


— e GAs 2) 9, (e242—99,G,) + e^ MG sin = MY nhn: (4.2) 
Via the substitution Y, > ec 3BAs—O)g the equation can be brought into schrédinger-like 


equation 
-Gr + Vain = MY Gn; (4.3) 

with the 5D effective schrodinger potential 
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4.2 Vector glueballs 
For vector gluebalsl Y, the 5D action is 
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The equation of motion for vector glueballs ¥ is 
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Via the substitution %⁄p > eas -0) y the equation can be brought into schrödinger-like 
equation 


-Hr + Vy hn = My Va, (4.7) 
with the 5D effective schrödinger potential 
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4.3 Tensor glueballs 
For tensor glueballs, the 5D action is 
1 j 
Sr =—5 f da/fge?( VihunVERMS — 2V Lh VN hym +2V mh NV yh 
-V mhV h + Mf (nv huy — h?)), (4.9) 

where h = gYNhmyn. With the constraint 

Vuh” =0, h=0, huw = eÆ, hmz=0, (4.10) 

The equation of motion for tensor glueballs “%,,, is 


— g B4A:-8) 9, (6841-99, 7) + eA M2, 9, = ep, Hr (4.11) 


Via the substitution #, > e 3 BAd) yy the equation can be brought into schrédinger- 
like equation 


-H + Vie Hy, = my Hr, (4.12) 
with the 5D effective schrodinger potential 
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4.4 Numerical results 


For numerical calculations, we have to fix parameters in the model. In the dynamical 
holographic model, there is only one free parameter, i.e., ug. We fix this parameter by 
fitting the scalar glueballs spectra from lattice results [2-5] as shown in Table 2. The 
lattice data in Table 2 indicates the slope of the Regge spectra is around 4GeV?, which is 
equivalent to ug ~ 1GeV in the dynamical holographic QCD model. 


n(0T*) Lat1 Lat2 Lat3 Lat4 Lat5 
No =3 Ne=3 Ne oo N: = Ne =3 
1 1475(30)(65)  1580(11) 1480(07)  1730(50)(80) — 1710(50)(80) 
2 2755(70)(120) © 2750(35) 2830(22) 2670(180)(130) 
3 3370(100)(150) 
4 3990(210)(180) 


Table 2. Lattice data for 0**glueball in unit of MeV. Lat1 data from Ref.[4], Lat2 and Lat3 data 
from Ref.[3], Lat4 [2] and Lat5 [5] are anisotropic results. 


We will also compare our results in the dynamical holographic QCD model with those 
in the hard-wall and soft-wall holographic QCD models. In the hard-wall holographic QCD 
model, the equation of motion for glueball @ is: 


o e489, (e°4° 0,G,) Zs eA: MG 5G, = MY nns (4.14) 


where c = 1 for vector glueballs and c = 3 for scalar and tensor glueballs. With UV 
boundary condition Y,,(€) = 0, the solution is: 


Gal) = 2 Jal Manz) (4.15) 


where n = ,/1+2c+c? + 4M4 5/2 and J is Bessel function. IR boundary condition 
0:Gn(Z2m) = 0 gives the discrete spectrum of the glueballs. Here zm is the hard cut-off, 
which is the only parameter in the hard-wall model, and can be fixed by the ground state of 
the scalar glueball 0++. When we take the mass for the lowest scalar glueball as 1730MeV, 
which fixes zm = 452MeV in the hard wall model. 


In the soft-wall model, with the dilaton background takes the quadratic form ® = ee 
and the metric structure is still AdSs, the Regge spectra for glueball can be derived as 
[24, 25] 


me = pd fante+14 (e+ ip rang, n=0,1,2, (4.16) 


where c = 3 in case of scalar and tensor and c = 1 in case of vector. There is also only one 


parameter ug in the soft-wall model. The Regge slope of 0** glueball gives ug = 1GeV 
(SW), and the the lowest scalar glueball mass 1730MeV gives ug = 0.5GeV (SWtt). 
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Table 3. The mass spectra of 0T* glueballs in the dynamical soft-wall model, compared with 
combined lattice data [2-5], and hard-wall model (zm = 452MeV), soft-wall model with SW} 
indicates ug is fixed by lowest mass of 0** glueball, and SWt indicates ug is fixed by Regge slope 
of OTF glueball. The unit is in MeV. 


Table 3 shows scalar glueball spectra in the dynamical soft-wall holographic QCD 
model, and the results are compared with combined lattice data [2-5] as well as hard-wall 
and soft-wall models, respectively. It is found that the hard-wall model cannot produce 
Regge spectra, and the soft-wall model cannot simultaneously produce the correct Regge 
slope and the ground state of the scalar glueball. As it was shown in Ref. [29], the Regge 
slope and the ground state of the scalar glueball can be correctly produced in the dynamical 
soft-wall holographic QCD model with only one parameter. 

With the parameter fixed by the scalar glueball, we calculate the vector and tensor 
glueball masses in the dynamical soft-wall model and compare with lattice data as well 
as results from hard-wall and soft-wall models. The results are shown in Table 4. It is 
observed that for vector and tensor glueballs, the results from the dynamical holographic 
QCD model are far away from lattice data, especially the masses for higher spin states 
are too heavy comparing with lattice data. The glueball massea from the hard-wall model 
are in general lighter than lattice results. Among the three models, the soft-wall model 
(SWtt) with the parameter fixed by the Regge slope can produce reasonable good results 
comparing with lattice data. However, all models cannot distinguish even and odd parity 
state for the glueball with the same spin. 

In the next section, we will improve the dynamical soft-wall holographic QCD model 
in order to produce reasonable glueball spectra. 
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Table 4. The mass spectra of vector and tensor glueballs in the dynamical soft-wall model, com- 
pared with lattice data, and hard-wall model (zm = 452MeV), soft-wall model with SW} indicates 
uG is fixed by lowest mass of Ott glueball, and SWẸł indicates ug is fixed by Regge slope of 0+ 
glueball. The unit is in MeV. 


5 Glueball spectra in modified dynamical soft-wall holographic QCD 
model 


As we observed from last section that the masses for higher spin glueballs are too heavy 
comparing with lattice data, while these states are reasonable in the soft-wall model. This 
indicates that only scalar glueballs are sensitive to the deformed metric, and other glueballs 
are not excited from this deformed metric background. Therefore we introduce a deformed 
5D mass squared for glueballs. In order to distinguish even and odd parity, we introduce 
the positive and negative coupling between the dilaton field and glueballs, respectively. 
With this set-up, now the 5D action for the scalar, vector and tensor glueballs 4 (x, z) take 
the following form: 
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Ss =- I dx /gse??( OuGING + MZ 5(2)2°), (5.1) 
1 1 

Sy = =r | Payme P SEMN + My 5(z)V”), (5.2) 
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where M2 (z) = Méee2#/8, p =1 for even parity and p = —1 for odd parity. 
The equation of motion for any glueball <& can be brought into schrédinger-like equa- 


tion 
— By + Va by, = my Dy, (5.4) 

with the 5D effective schrodinger potential 
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where c = 1 for 1-form and c = 3 for 0-form and 2-form, and M2, z is the value given in 
Table 1. 
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Figure 1. The effective schrodinger potential V of scalar and vector glueballs in the soft-wall model 
(green thick line), dynamical soft-wall model (orange dashed line) and dynamical soft-wall model 
with modified M?(blue line). 
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Figure 2. The effective schrédinger potential V of tensor glueballs in the soft-wall model (green 
thick line), dynamical soft-wall model (orange dashed line) and dynamical soft-wall model with 
modified M? (blue line). 


Compare the 5D effective schrédinger potential Eq.(5.5) with Eqs.(4.4), (4.8) and 
(4.13), we can see the effect of the deformed 5D mass square M?(z) = M?e~2®/8 is to 
counteract the deform metric background. In Fig. 1 and Fig. 2, we show the 5D effec- 
tive schrédinger potential Eq.(5.5) as a function of z and compare with results from the 
soft-wall model and the original dynamical soft-wall model. It is found that at infrared 
(IR), except the scalar glueball 0**, the 5D effective schrédinger potential for other glue- 
balls in the modified dynamical soft-wall holographic QCD model coincide with those from 


soft-wall model. The parity difference p = + only brings the difference of the 5D effective 
schrédinger potential in the range of 0.5 < z < 2. 

The final results of the glueball spectra in the modified dynamical holographic QCD 
model are shown in Fig. 3 and in Table 5 with details. It is found that with only one 
parameter ug = 1GeV, which is fixed by the Regge slope of the scalar gluball spectra, 
one can produce other glueballs spectra agree well with lattice data, except three trigluon 
glueball states 0~~, 0+7 and 2+7, whose masses are 1.5 GeV lighter than lattice results. 


=] i= 


Considering that we only take the simplest quadratic dilaton profile, which corresponds 
to dimension-2 gluon condensate (or effectively two-gluon condensate) in the vacuum, our 
results might indicate that these three trigluon glueballs 077, 0+7 and 2*~ are dominated 
by three-gluon condensate contribution. 
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Figure 3. The mass spectra of glueballs in the modified dynamical soft-wall model with u = 1 GeV 
(line) compared with the lattice result [2-6] (rectangle). 


6 Conclusion and discussion 


In this work, we study scalar, vector and tensor glueball spectra in the framework of 5- 
dimension dynamical holographic QCD model, where the metric structure is deformed 
self-consistently by the dilaton field. It is found that only scalar glueballs are excited 
from this deformed metric background, and other glueballs excited from this deformed 
metric background are much heavy comparing with lattice data. Therefore, for higher spin 
glueballs, we introduce a deformed 5-dimension mass in order to counteract the effect of the 
deformed metric background. In order to distinguish glueballs with even and odd parities, 
we introduce the positive and negative coupling between the dilaton field and glueballs. 
With these set-ups, we calculate the glueball spectra with only one free parameter in 
the dynamical holographic QCD model, which is fixed by the scalar glueball spectra. It 
is found that all two-gluon glueball spectra produced in the dynamical holographic QCD 
model are in good agreement with lattice data. We investigate six trigluon glueballs, among 


these trigluon glueballs, the produced masses for 1*~ and 277 are in good agreement with 
lattice data, and the produced masses for 077, 0+7 and 2*7 are around 1.5 GeV lighter 
than lattice results. Considering that we only take the simplest quadratic dilaton profile, 
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Table 5. The mass of glueball spectra in Lattice QCD [2-6], Flux tube model [7], QCDSR [9, 10, 
35, 37-39] and modified dynamical soft-wall model. Note that 0T§ is trigluonium. The unit is in 
GeV. 


which corresponds to dimension-2 gluon condensate (or effectively two-gluon condensate) 
in the vacuum, our results might indicate that the three trighion glueballs 077, 0*~ and 
2*~ are dominated by three-gluon condensate contribution. Further studies with more 
complicated dilaton profile are needed. 
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